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Abstract 



A proposal to test the proposed time-reversal and inversion breaking phase 
in the Pseudogap region of the Cuprate compounds through the variation of 
Second-harmonic generation intensity with temperature and polarization and 
angle of incidence is presented. 
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I. INTRODUCTION 



Feibig et al. 1 have discussed the conditions for second harmonic generation proportional 
to the Time-reversal breaking order parameter for the case of magnetoelectric symmetry 2 , 
i.e. for the case that Time-reversal Symmetry R is broken and so is Inversion I but their 
product RI is preserved. Following earlier suggestions 3 ' 4 , ARPES experiments with circularly 
polarized photons 5 have detected such a phase in the so-called pseudogap phase in the phase 
diagram of a high temperature superconducting compound of the Copper-Oxide (Cu-O) 
family. Time-reversal is broken due to circulating currents in specific patterns within the 
unit cell with translational symmetry preserved. If this discovery is correct, it togather with 
the extensive evidence for Quantum Critical properties in the normal phase 6 around the 
doping- density for the highest T c implies that the basic framework for the theory of the 
Cuprates is in place. It is therefore very important to verify the discovery by independent 
experiments. Accordingly, we work out here the symmetry considerations for magneto- 
electric induced second harmonic generation in the Cuprates. 

Many terms in the multipole susceptibility contributes to the second harmonic gener- 
ation (SHG). The relevant contributions for the processes discussed here come from the 
magnetic-dipole, Mk(2u) oc Xijk( u i cu, 2uj)Ei{uj)Ej{uj), and the electric-dipole, Ek(2u) oc 
Xtjk( u i °°i 2uj)Ei(uj)Ej(uj). While the former is generally present in the symmetry of the 
cuprate compounds, the latter occurs in the predicted pseudogap phase with xtjk transform- 
ing as an odd tensor under time-reversal. The intereference effects due to the simultaneous 
presence of the two-time-reversal odd tensors xTjk an d Xtjk are responsible for the effects 
discussed here. 



II. SYMMETRY OF THE NORMAL PHASE 

Most Cu — O compounds possess a centre of symmetry in the unit-cell. We consider 
the body-centered tetragonal crystals with symmetry group 4/mmm(D 4 / l ), denoting the 
presence of a four-fold axis and a mirror plane orthogonal to this axis and two other mirror 
planes or orthorhombic crystals with symmetry group mmm (D2h), (three orthogonal mirror 
planes only). Due to the centre of inversion, the third rank polar tensor xtjk arising from 
electric- dipole effects is zero but the third rank axial tensor xTjk arising from magnetic-dipole 
effects is non-zero. In other words a term in the Free-energy $ of the form 

$ = xT jk ErE 3 H k (1) 

exists. Correspondingly a magnetization for applied electric fields 

M k = x^E.Ej (2) 

is non-zero. 

We specify the tensors in this paper using the following co-ordinates: The c axis of 
the crystal coincides with the z direction, a and b crystalline axis (those along the nearest 
neighbor Cu-Cu vector) are along the (x+y) / y/2 and (x — y) / ^J2 directions. For the A/mmm 
group, three independent components of xTjk are non-zero 7 . The vector M takes the form: 
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M 



X?E Z E X 




(3) 



X? = UZy + XZy) (4) 



For the mmm group, six independent elements of xTjk ex ist and all three cartesian 
components of M are bilinear in E's. We present detailed results for the 4/ 'mmm crystalline 
symmetry. Very similar conclusions hold for the mmm crystalline symmetry. 



III. TIME-REVERSAL BREAKING STATE 

Two different phases breaking time-reversal are shown to be possible in the pseudogap 
phase ( 3 , 4 ). The current patterns in these phases in the Cu-0 planes are shown in Fig. (1). 

One called the 0// phase has the symmetry mmm (Fig. 1 b). This group is in the 
Magnetoelectric class, where both time-reversal (R) and inversion (7) are broken but their 
product is preserved. The other possible phase 0/ (Fig. 1 a) preserves inversion while 
breaking Time- reversal. It has the symmetry 4/ mmm. So the effects discussed here do not 
occur in such a phase. 

As will be clear from Fig. (lb), the point group mmm has the following elements: the 
identity, two-fold axis rotation around y (C y ), the reflections in the x — y and y — z planes 
(<t 2 and a x ), RI, Ra y , RC X and RC Z . 

In the group mmm, axial tensors of third rank which are invariant under time-reversal 
Xijk^ are allowed 7 . The independent components of M are 



X™E y E z 

M = I X yE x E z I . (5) 



XTE X E 



v 



where xjki + Xkji = XT an< ^ k run over x > Vi z - But more importantly, polar-tensors of 
odd-rank xtjk > which change sign under time- reversal are also allowed. So a Polarization 
vector of the form 



P = 



( xf ] E x E y \ 

X^El + x^El + x^X 
\ xf ] E z E y J 



(6) 



is allowed. Here xt^ — Xyxx + Xxyl and xt^ = Xyzz + Xzyz- These x' s are expected to be 
proportional to the order parameter Q n . 

From these we can construct the relevant part of the Poynting vector S: 

d 2 P „ dM 

S = ^ + VX ^ (7) 
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FIGURES 
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FIG. 1. Current patterns for the time-reversal breaking states proposed in Ref. 3 ' 4 . 

The SHG intensity is I oc \S\ 2 . As shown below, interesting effects arise in experimental 
geometries in which both an d X m are used and interfere. Then effects linear in the 
magnetoelectric order parameter are observable. If light is propogated in the sample along 
any of the axes, V x M = 0. Then only the electric-dipole susceptibility contributes to the 
SHG and I oc \x e \ 2 which is in turn proportional only to the square of the order parameter. 

We fix the axes with respect to the crystal as in Fig. (2). For an arbitrary direction 
of incidence, Figure 2 shows the triad form by the propagation vector q (specified by the 
polar angle 9 and the azimuthal angle 0) , and the components of the polarization along the 
orthogonal directions e xy in the x-y plane and e qz in the q-z plane. To maximize the second 
term in Eq. (7) we chose a polar incidence angle 9 = tt/4. 
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FIG. 2. Geometry of incident light vector q. The axes are fixed to the crystal: x and y are 
rotated 7r/4 with respect to the crystalline axis a and b as explained in the text. 



For light linearly polarized along the direction E = Ee 



qz: 



oj sin 0s 



q|(cos 2 <^ 



E oj I u) cos 0s% — |q| sin y, , Vr , v „ 

cusm^st — Iqlfcos 2 (fry™ + sin 2 0x™)/2 



(8) 



where, 



si 
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(cos0(xi 



e(c) 



e(c) 

X 2 



cos 2 0xS - sin 2 x£ - xS)/(2V2) 



/ / / e(c) e(c) 
(COS0(X1 + X 2 



xxi/ " t" *\.yyy <\.zzy , 
_|_ air.2 ^ v e(c) _ e(c). 



,, ( + Bin^(C-X e r ; ))/2 



cos 2 0xS - sin 2 xS - xS)/(2V2) 



(9) 



The component are given in the basis (q, e xy , e qZ/ 



I oc E 4 u 2 [cj 2 (sin 2 0(|s^| 2 + |^| 2 ) + cos 2 0| s 



if) 



H\\ sin0cos0( X ™ + Xy + XDI78 + I cos' 
u\q\ sin0(cos 2 + Xy + xTY 



iXy + sin 2 . 
3?Sg (cos 2 . 



,m\2\ 



(10) 



</) X r + sin 2 X ™)*)/2 



The x's denned below are complex and x* denotes complex conjugation. The last term in 
Eq. (10) is linear in both \ m and x e an d is therefore proportional to the order parameter. 
For an incident light linearly polarized, E = Ee xy , 



oj sin <ps\ 

E £ oj | oj cos (j)S2 + | q| sin (j) cos 4>x™ / V% 



(11) 



oj sm 0s\ 



where 
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[sm cos <; 

,2 A{~A C ) 



etc) ■ 2 

i ' — sm i 



-(sin 2 0( X r ; + X:g) + cos 2 



Ac) 
^yyy 



Ac) 

^.xxy 



(12) 
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It is interesting to consider circularly polarized light with right E + = —l/y/2(E xy + iE zq ) 
and left E_ = l/y/2(E xy — iE zq ) polarizations and evaluate the intensity / oc |S| 2 for each 
case. This gives 

/ oc \E ± \ V (o;V(0) e | 2 + |q| V(0)"T ± cos(0)i?e^|q|^(0) e (^ni(0) m )*)) (13) 

Again the last term is the interesting term as it is linear in both \ m an d X e an d is therefore 
proportional to the order parameter. Its sign is changed on reversing the polarization of 
the light or by reversing time (i.e. by reversing the domains or reversing the direction of 
propagation for a fixed domain i.e. changing <fi to n + <fi). For cos0 = ±1, that is when q 
belongs to the broken symmetry plane y, the difference between both circular polarization 
is maximized. For this case, 

1^(0)1 

2\m m \ 

^(7T/2) e ( Xmt (7T/2) m ) 

In summary we have determined the SHG intensities, as a function of polarization and 
direccion of the incident light, for different group symmetries corresponding to the normal 
and time-reversal breaking symmetry states. In all the states there is a contribution to the 
SHG coming from the magnetic-dipole terms but the electric-dipole contribution appears 
only when the inversion and time-reversal are broken with their product preserved, as in the 
Qn phase. The interference of these two terms produce linear terms in both x m and x e an d 
therefore proportional to the order parameter. This term change sign reversing the domains 
On to —On or changing the direction of incident light by the operation q — > o x o y <\ (which 
has the effect of changing^ to <fi + ir). The interference term also changes sign on changing 
from left to right circularly polarized light. 

The presence of many different 0// domains within the spot size of the incident beam 
will blur this effect. ARPES experiments 5 indicate domains to be of 0(100 microns). 

The proposed experiment are a rather thorough test of the symmetries of the predicted 
phase because they rely on broken R and a y (or cr x )) but their product preserved below 
the pseudogap temperature in the cuprates. Breaking reflection symmetry alone would not 
produce the predicted effects in Second-harmonic generation, while the ARPES experiment 5 
is in principle compatible 4 with breaking some reflection symmetries alone. Some other 
experiments using resonant x-ray absorption and diffraction have also been proposed ( 8 ). 

Acknowledgments: We wish to thank Professor Harry Tom for discussion of the feasi- 
bility of experiments and Professor J. Orenstein who asked the questions which led to this 
investigation. 



|X2 I IXl I \Xyyy I ^ \Xxxy Xzzy I (14) 

= | Y m | 2 + | Y m + v m | 2 
I Ay I 1 \A.x 1 A.z I 

= (xS + xS - 2x£)(xr + x™y + (xi (c) + X3 (c) )(x2-r 
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